A new integral representation for quasi-periodic scattering problems in two dimensions

نویسندگان

  • Alex Barnett
  • Leslie Greengard
چکیده

Boundary integral equations are an important class of methods for acoustic and electromagnetic scattering from periodic arrays of obstacles. For piecewise homogeneous materials, they discretize the interface alone and can achieve high order accuracy in complicated geometries. They also satisfy the radiation condition for the scattered field, avoiding the need for artificial boundary conditions on a truncated computational domain. By using the quasi-periodic Green’s function, appropriate boundary conditions are automatically satisfied on the boundary of the unit cell. There are two drawbacks to this approach: i) the quasi-periodic Green’s function diverges for parameter families known as Wood’s anomalies, even though the scattering problem remains well-posed, and ii) the lattice sum representation of the quasi-periodic Green’s function converges in a disc, becoming unwieldy when obstacles have high aspect ratio. In this paper, we bypass both problems by means of a new integral representation that relies on the free-space Green’s function alone, adding auxiliary layer potentials on the boundary of the unit cell strip while expanding the linear system to enforce quasi-periodicity. Summing nearby images directly leaves auxiliary densities that are smooth, hence easily represented in the Fourier domain using Sommerfeld integrals. Wood’s anomalies are handled analytically by deformation of the Sommerfeld contour. The resulting integral equation is of the second kind and achieves spectral accuracy. Because of our image structure, inclusions which intersect the unit cell walls are handled easily and automatically. We include an implementation and simple code example with a freely-available MATLAB toolbox. A. H. Barnett Department of Mathematics, Dartmouth College, Hanover, NH, 03755, USA. tel:+1-603646-3178. fax:+1-603-646-1312 E-mail: [email protected] L. Greengard Courant Institute, New York University, 251 Mercer St, NY, 10012, USA. E-mail: [email protected] 2 Alex Barnett, Leslie Greengard

برای دانلود رایگان متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Challenges in Geometry , Analysis and Computation : High Dimensional Synthesis

We present a new technique, based on a clever numerical minimization of certain Sobolev norms, for the numerical solution of partial differential equations. The new method is very high order even for complex geometries and variable coefficient problems. One surprising feature of the method is that we prefer to pose the problem in first-order form. The system does not have to be square as long a...

متن کامل

Robust high-order numerical scattering from multi-layer dielectric gratings using a new integral representation for quasi-periodic fields

Many numerical problems arising in modern photonic and electromagnetic applications involve the scattering of a plane wave from a piecewise-homogeneous medium, periodic in one direction. Boundary integral equations are an efficient approach to solving such boundary-value problems with high-order convergence. The standard way to periodize is then to replace the free-space Green’s function kernel...

متن کامل

Efficient Solution of Acoustic and Electromagnetic Scattering Problems in Three-Dimensional Periodic Media

We present an accurate and efficient numerical method, based on integral Nyström discretizations, for the solution of three dimensional wave propagation problems in piece-wise homogeneous media that have two-dimensional (in-plane) periodicity (e.g. photonic crystal slabs). Our approach uses (1) A fast, high-order algorithm for evaluation of singular integral operators on surfaces in three-dimen...

متن کامل

Robust fast direct integral equation solver for quasi-periodic scattering problems with a large number of layers.

We present a new boundary integral formulation for time-harmonic wave diffraction from two-dimensional structures with many layers of arbitrary periodic shape, such as multilayer dielectric gratings in TM polarization. Our scheme is robust at all scattering parameters, unlike the conventional quasi-periodic Green's function method which fails whenever any of the layers approaches a Wood anomaly...

متن کامل

A fast direct solver for quasi-periodic scattering problems

We consider the numerical solution of the scattering of time-harmonic plane waves from an infinite periodic array of reflection or transmission obstacles in a homogeneous background medium, in two dimensions. Boundary integral formulations are ideal since they reduce the problem to N unknowns on the obstacle boundary. However, for complex geometries and/or higher frequencies the resulting dense...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

عنوان ژورنال:

دوره   شماره 

صفحات  -

تاریخ انتشار 2010